Anisotropy at the end of the cosmic ray spectrum? 
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The starburst galaxies M82 and NGC253 have been proposed as the primary sources of cosmic 
rays with energies above 10 18 7 eV. For energies > 10 20 ' 3 eV the model predicts strong anisotropies. 
We calculate the probabilities that the latter can be due to chance occurrence. For the highest 
energy cosmic ray events in this energy region, we find that the observed directionality has less 
than 1% probability of occurring due to random fluctuations. Moreover, during the first 5 years of 
operation at Auger, the observation of even half the predicted anisotropy has a probability of less 
than 10~ to occur by chance fluctuation. Thus, this model can be subject to test at very small 
cost to the Auger priors budget and, whatever the outcome of that test, valuable information on 
the Galactic magnetic field will be obtained. 



INTRODUCTION 

Soon after the microwave echo of the big bang was 
discovered, Greisen, Zatsepin, and Kuzmin (GZK) noted 
that the relic photons make the universe opaque to cos- 
mic rays (CRs) of sufficiently high energy 0. This oc- 
curs, for instance, for protons with energies beyond the 
photopion production threshold (A(1232) resonance). 
After pion production, the proton (or perhaps, instead, 
a neutron) emerges with at least 50% of the incom- 
ing energy. A similar phenomenon (of energy degra- 
dation) occurs for nuclei due to processes of photodis- 
integration. Therefore, the characteristic attenuation 
length for extremely high energy (10 20 eV < E < 
10 20 ' 5 eV) hadrons is less than 100 Mpc, decreasing down 
to 10 Mpc with rising energy 0]. The survival proba- 
bility for extremely high energy (EHE) 7-rays (propa- 
gating on magnetic fields 3> 10~ n G) to a distance d, 
P(> d) ~ exp[— d/6.6 Mpc], becomes less than 10~ 4 af- 
ter traversing a distance of 50 Mpc . This implies that 
the GZK sphere -4] represents a small fraction of the size 
of the universe. Consequently, if the CR sources are uni- 
versal in origin, the energy spectrum should not extend 
(except at greatly reduced intensity) beyond ~ 10 20 eV, a 
phenomenon! known as the GZK cutoff. Even though the 
Haverah Park 0, Yakutsk |(J, Fly's Eye 0, and HiRes Q 
data show statistically significant evidence for such a cut- 
off (more than 5cr independent of the sample used as 
a basis for extrapolation), the AG AS A ground array de- 
tected a handful of events with energies > 10 20 eV [loj . 
as opposed to about 2 expected from the GZK cutoff. 
Moreover, within statistical uncertainty (which is large 
above 10 20 eV) the flux of CRs above 10 187 eV reported 
by the AGASA Collaboration [Jjj is consistent with a 
E~ 2 ' 7 spectrum up to the highest observed energies, sug- 
gesting that a single acceleration mechanism is responsi- 
ble for all the events beyond that energy, unless of course 
a very unlikely matching of spectra can account for the 
smoothness of the CR energy distribution. 

In order to analyze the effect of energy losses on the 
observed spectrum, it is convenient to introduce the ac- 



cumulation factor / acc , defined as the ratio of energy- 
weighted fluxes for "low" (10 187 eV - 10 19 5 eV) and 
EHECRs. With this in mind, if the Earth is located in 
a typical environment and all CR-sources have smooth 
emission spectra, the observed spectrum above 10 187 eV 
should have an offset in normalization between low and 
EHE given by / acc . For CR protons and nuclei with 
uniform distribution of sources active over cosmological 
times, the cutoff due to photopion and photodisintegra- 
tion processes relates the accumulation factor to a ra- 
tio of attenuation lengths and leads to f arr ~ 100. 
The smoothness of the observed CR spectrum [lTj, viz. 
/acc ~ 1, seems to indicate that the power of nearby 
sources must be comparable to that of all other sources 
(redshift z > 0.5) added together. 

The simplest explanation, i.e., nearby sources are sig- 
nificantly more concentrated, does not seem to be the 
case. Specifically, if one simply assumes that the dis- 
tribution of CR sources follows the distribution of nor- 
mal galaxies, the local overdensity is only a factor of two 
above the mean, and thus insufficient to explain the mea- 
sured flux above 10 20 eV Furthermore, the arrival 
direction of the super-GZK events is consistent with an 
isotropic distribution of sources (even when some level of 
clustering was already detected |ljj), in sharp contrast to 
the anisotropic distribution of light within 100 Mpc [lij . 
A way to avoid the problems with finding plausible as- 
trophysical explanations is to look for solutions involving 
physics beyond the standard model . While the invo- 
cation of such new physics is an intringuing idea, there 
are now constraints that call into question the plausibility 
of some of these ideas [lrij . 

Recently, it was suggested that the observed near- 
isotropy of arrival directions could be due to a diffuse 
propagation of EHECRs I n this work, we exam- 

ine specific candidate sources for this hypothesis. These 
are the starburst galaxies M82 and NGC253 which have 
been shown to reproduce the main features of the ob- 
served flux I n particular, we study here the critical 
aspect of a residual anisotropy that emerge beyond the 
GZK energy limit after deflection in Galactic and ex- 
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tragalactic magnetic fields. Specifically, we estimate the 
probability that an apparent correlation between the ar- 
rival directions of the highest energy events and the two 
starbursts can originate as a purely random fluctuation. 
After that, we study the sensitivity of Auger Observatory 
to the model. 



DIFFUSE PROPAGATION OF COSMIC RAYS IN 
A MAGNETIZED NEIGHBORHOOD OF THE 
GALAXY 

A popular explanation considered recently |l7| for an 
isotropic distribution of arrival directions entertains the 
existence of large scale intervening magnetic fields, so 
that even EHECRs propagate diffusively. Indeed, there 
are some measurements of diffuse radio emission from 
the bridge area between Coma and Abell superclusters 
that under assumptions of equipartition allows an esti- 
mate of 0.2 — 0.6 /iG for the magnetic field in this re- 
gion [J^. Such a strong magnetic field (which is com- 
patible with existing upper limits on Faraday rotation 
measurements 20]) could be possibly understood if the 
bridge region lies along a filament or sheet of large scale 
structures [U- In light of this, it appears plausible, 
though subject to verification, to assume that our Local 
Supercluster contains a large scale magnetic field (say, 
10" 8 G < B < IQ- 6 G which provides sufficient 

bending to EHECR orbits, camouflaging the exact loca- 
tion of the sources. 

Diffusion has two distinctive regimes. Particles that 
are trapped inside magnetic subdomains (of size £-wipc = 
£/Mpc) follow Kolmogorov diffusion. In such a case, the 
functional dependence of energy of the difussion coeffi- 
cient (for protons) is found to be [23| 
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where B^q is the magnetic field strength in units of /iG 
and E20 is the particle's energy in units of 10 20 eV. With 
rising energy, the Larmor radius of the particles starts 
approaching £ and there is a transition to Bohm diffusion. 
The diffusion coefficient in this regime is of order the 
Larmor radius times velocity (~ c). 

If CRs propagate diffusively, the radius of the sphere 
for potential proton sources becomes significantly re- 
duced. This is because one expects negligible contribu- 
tion to the flux from times prior to the arrival time of 
the diffusion front, and so the average time delay in the 
low energy region, 

d 2 

Tdclay « , (2) 

must be smaller than the age of the source, or else the 
age of the universe (if no source within the GZK ra- 
dius is active today, but such sources have been active 



in the past). Note that the diffuse propagation of EHE 
protons requires magnetic fields ~ 1/xG. Therefore, for 
typical coherence lengths of extragalactic magnetic fields 
{£ ~ 1 Mpc) the time delay of CRs with E « 10 187 eV 
cannot exceed Td c i ay ^ 14 Gyr, yielding a radius of 
d ~ 30 Mpc. In the case CR sources are active today, the 
radius for potential sources is even smaller d ~ 5 Mpc. 
Centaurus A, at a distance of 3.4 Mpc and galactic co- 
ordinates I = 310°, b = 20°, is the nearest active galaxy, 
and the only one within a distance of 5 Mpc. Phcnomeno- 
logical arguments identify Centaurus A as a plausible 
progenitor of all CRs observed on Earth with energies 
> 10 18 - 7 eV However, detailed numerical simu- 

lations seem to indicate that large scale magnetic fields 
0(/iG) cannot provide sufficient angular deflection to ex- 
plain all the observational data: (1) the large deflection 
angle of the highest energy event recorded by the Fly's 
Eye experiment (see Table I) with respect to the line of 
sight to Centaurus A must be explained as a 2a fluc- 
tuation jl^l, (2) for an emission spectrum cx E~ 2A and 
maximum injection energy of 10 21 eV, the angular power 
spectrum shows a 3ct quadrupole deviation from AGASA 
observations [30| . 

If magnetic fields in the nanogauss range exist in the 
neighborhood of the Galaxy, it is possible that ultrahigh 
energy cosmic ray nuclei could diffuse sufficiently in order 
to attain the observed near-isotropy. For a CR nucleus 
of charge Ze in a magnetic field B n Q = B/10~ 9 G, the 
Larmor radius is 
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In this case, the sphere of potential sources is severely 
constrained by the GZK cutoff: less than 1% of iron nu- 
clei (or any surviving fragment of their spallations) can 
survive more than 3 x 10 14 s with an energy > 10 20 ' 5 eV. 
Therefore, the assumption that EHECRs are heavy nu- 
clei implies ordered extragalactic magnetic fields B n Q < 
15 — 20, or else nuclei would be trapped inside magnetic 
subdomains suffering catastrophic spallations. There are 
two candidate sources within the GZK sphere; namely, 
the nearby (d ~ 3 Mpc) metally-rich starburst galax- 
ies M82 (I = 141°, b = 41°) and NGC253 (I = 89°, 
b = —88°) [3l|| . Phenomenological considerations based 
on analytical estimates of the diffusion coefficient and ap- 
proximations to the photodisintegration losses and angu- 
lar deflections suggest that the power of these starbursts 
is enough to provide all CRs observed on Earth above 
10 18 ' 7 eV (lif. This analytical study is consistent with 
Monte Carlo simulations |33|. Specifically, the spectrum 
observed by AGASA can be fitted with a single source lo- 
cated at d = 3.2 Mpc if the spectrum of nuclei is oc E^ 1 - 6 . 
This is a hard spectrum compared to the expected E~ 2 
from the Fermi mechanism. In this context, it should be 
noted that the fit in has strong statistical weight from 
points near 10 19 eV. However, there is significant system- 
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TABLE I: The highest energy cosmic rays. The energy reso- 
lution for the AGASA experiment was taken from Ref. 

El- 
Date Experiment E[EcV] I b Ref. 
89/05/07 Yakutsk 300t\ 7 s 162.0° 2.0° [25] 
91/10/15 Fly's Eye 320 ± 90 163.4° 9.6° [261 
93/12/03 AGASA 213 ± 75 130.5° -41.4° [27] 
01/05/10 AGASA 280 ± 98 106.3° -39.0° [28] 



atic uncertainty in the observed energy spectrum in this 
region. For example, as recently noted a downward 
shift of 11% in the AGASA energy calibration is required 
in order to bring the resulting spectrum into agreement 
with Fly's Eye data. This softening of the observed ener- 
gies will require a steeper Fermi-like injection spectrum. 

The most salient feature of the starburst hypothesis is 
the prediction of an anisotropy at the high end of the 
spectrum. For an extragalactic, smooth, magnetic field 
of w 15 — 20 nG, diffusive propagation of particles be- 
low 10 20 eV evolves to nearly complete isotropy in the 
CR arrival directions fl8t l32| . Above this critical energy 
there is a transition range (up to 10 20 3 eV) where the 
combined bending in extragalactic and Galactic magnetic 
fields leads to loss of directionality. With rising energy, 
the average deflection in the extragalactic magnetic field 
is significantly reduced, and is roughly 10° < 9 < 20° 
in the energy range 10 20 3 eV < E < 10 20 5 eV 0. In 
order to incorporate typical uncertainties in energy reso- 
lution as well as those in the Monte Carlo simulation , 
we will increase the upper limit of this deflection to 30°. 
Heavy nuclei suffer additional deflection in the Galactic 
magnetic field. 

The large scale structure of the Galactic magnetic field 
carries substantial uncertainties, because the position of 
the solar system does not allow global measurements. 
The average field strength can be directly determined 
from pulsar observations of the rotation and dispersion 
measures average along the line of sight to the pulsar 
with a weight proportional to the local free electron den- 
sity, (Bn) w 2^tG |3^| . Measurements of polarized syn- 
chrotron radiation as well as Faraday rotation of the 
radiation emitted from pulsars and extragalactic radio 
sources revealed that the global structure of the magnetic 
field in the disk of our Galaxy could be well described by 
spiral fields with 27r (axisymmetric, ASS) or tt (bisym- 
metric, BSS) symmetry |34j. In the direction perpen- 
dicular to the Galactic plane the fields are either of odd 
(dipole type) or even (quadrupole type) parity. Discrimi- 
nation between these models is complicated. Field rever- 
sals are certainly observed (in the Crux-Scutum arm at 
5.5 kpc from the Galactic center, the Carina-Sagittarius 
arm at 6.5 kpc, the Perseus arm at 10 kpc, and pos- 
sibly another beyond H^)- However, as discussed by 



Vallee 36], turbulent dynamo theory can explain field 
reversals at distances up to ~ 15 kpc within the ASS 
configuration. Interestingly, if the Galactic field is of the 
ASS type, CRs entering the Galaxy with I < 180° are 
deflected towards increasing values of I and decreasing 
values of \b\ H3- Consequently, as we show in what fol- 
lows, each arrival direction given in Table I can be traced 
backwards to one of the starbursts. 

The field strength in the Galactic plane [z = 0) for the 
ASS model is generally described by [3?| 

B{p,6) = B (p) cos 2 [9- PMp/Zo)} , (4) 

where 9 is the azimuthal coordinate around the Galactic 
center (clockwise as seen from the north Galactic pole), 
p is the galactocentric radial cylindrical coordinate, and 

B (P) = — tanh 3 (p/ Pl ) pG . (5) 
P 

Here, £0 = 10.55 kpc stands for the galactocentric dis- 
tance of the maximum of the field in our spiral arm, f3 = 
1/ tanp (with the pitch angle, p = —10°), r — 8.5 kpc is 
the Sun's distance to the Galactic center, and p\ = 2 kpc. 
The 9 and p coordinates of the field are correspondingly, 

B 6 = B(p, 9) cos p , B p = B(p, 9) ship . (6) 

The field strength above and below the Galactic plane 
(i.e., the dependence on z) has a contribution coming 
from the disk and another from the halo, 

B(p,9,z) = B(p,9) tanh(z/z 3 ) 



2cosh(z/zi) 2cosh(z/z2) 




where z\ = 0.3 kpc, zi = 4 kpc and Z3 = 20 pc. Fig- 
ure n shows the extent to which the observed arrival di- 
rections of the CRs listed in Table I deviate from their 
incoming directions at the Galactic halo because of bend- 
ing in the magnetic field given in Eq. J7J). The incom- 
ing CR trajectories are traced backwards up to distances 
of 20 kpc away from the Galactic center, where the ef- 
fects of the magnetic field is negligible. The diamond at 
the head of each solid line denotes the observed arrival 
direction, and the points along these lines indicate the 
direction from which different nuclear species (with in- 
creasing mass) entered the Galactic halo. In particular, 
the tip of the arrows correspond to incoming directions 
at the halo for iron nuclei, whereas the circles correspond 
to nuclei of neon. Additionally shown in the figure, in- 
dicated by stars, is the location of the two starbursts. 
Regions within the dashed lines comprise directions ly- 
ing within 20° and 30° degrees of the starbursts. It is 
seen that trajectories for CR nuclei with Z > 10 can be 
further traced back to one of the starbursts, within the 
uncertainty of the extragalactic deviation. 
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FIG. 1: Directions in Galactic coordinates of the four highest 
energy cosmic rays at the boundary of the Galactic halo. The 
diamonds represent the observed incoming directions. The 
circles and arrows show the directions of neon and iron nuclei, 
respectively, before deflection by the Galactic magnetic field. 
The solid line is the locus of incoming directions at the halo 
for other species with intermediate atomic number. The stars 
denote the positions of M82 and NGC253. The dashed lines 
are projections in the (l,b) coordinates of angular directions 
within 20° and 30° of the starbursts. 



The trajectories in Fig. ^ result from motion in the 
regular component of the Galactic magnetic field. How- 
ever, there is some evidence supporting the existence of 
a random component roughly comparable in magnitude 
to the regular component pq . Thus, using the random 
walk formulation [39j with coherence lengths of ~ 1 kpc, 
we estimate that the trajectories should be broadened by 
an angle Si3 ~ 40° (Z/20). The effect of this broadening 
on our analysis will be further discussed in the following 
section. 



The effects of the BSS configuration are completely dif- 
ferent. Because of the averaging over the frequent field 
reversals, the resulting deviations of the CR trajectories 
are markedly smaller, and in the wrong direction for cor- 
relation of current data with the starburst sources. We 
note that the energy-ordered 2D correlation distribution 
of the AGASA data is in disagreement with expectations 
for positively charged particles and the BSS configura- 
tion m. 



ASSESSMENT OF RANDOM COINCIDENCES 

We now attempt to assess to what extent these corre- 
lations are consistent with chance coincidence. To do so, 
we first observe that the angular deviation of a CR arrival 
direction at the outer edge of the galaxy with respect to 
the straight line of sight is roughly 
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Hence, for rf k 5—10 kpc, the average deflection of 

heavy nuclei with energies in the range 10 20 3 eV < E < 
10 20.5 eV jg 30 o < £ < 4[r We arrive at the effective 

angular size of the source in a two-step process. Before 
correcting for bias due to the coherent structure of the 
Galactic magnetic field, the deflections in the extragalac- 
tic and Galactic fields (regular and random components) 
may be assumed to add in quadrature, so that the an- 
gular sizes of the two sources are initially taken as cones 
with opening half- angles between 40° and 60°, which for 
the purpose of our numerical estimate we approximate 
to 50°. However, the global structure of the field will 
introduce a strong bias in the cosmic ray trajectories, 
substantially diminishing the effective solid angle. The 
combined deflections in the I and b coordinates mentioned 
above concentrate the effective angular size of the source 
to a considerably smaller solid angle j3?| ■ As a conserva- 
tive estimate, we retain 25% of this cone as the effective 
source size. A clear prediction of this consideration is 
that the incoming flux shows a strong dipole anisotropy 
in the harmonic decomposition. 

In order to assess the likelihood of a random occurrence 
of the predicted spatial distribution we performed numer- 
ical simulations in the spirit of Ref. 41.]. By randomly 
generating four CR positions in the portion of the sky 
accessible to the existing experiments (declination range 
S > —10°), an expected number of random coincidences 
can be obtained. The term "coincidence" is herein used 
to label a synthetic CR whose position in the sky lies 
within an effective solid angle fl c s of either starburst. 
f2 e ff is characterized by a cone with opening half-angle 
reduced from 50° to 24° to account for the 75% reduc- 
tion in effective source size due to the magnetic biasing 
discussed above. Cosmic ray positional errors were con- 
sidered as circles of 1.6° radius for AGASA. For the other 
experiments the asymmetric directional uncertainty was 
represented by a circle with radius equal to the average 
experimental error. Figure [3 presents the simulation re- 
sults. There are, as we have seen, four real coincidences. 
However, the random prediction for the mean number of 
coincidences is 0.81 ±0.01 423 . The Poisson probabil- 
ity 0| for the real result to be no more than the tail of 
the random distribution is just 
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FIG. 2: Simulation results for 4500 trials (a larger number 
of trials do not modify these results). 99% of the simula- 
tions yield results with less than 4 random coincidences, as 
described in the text. 



where A is the mean value of the random results. Al- 
ternatively, we may analyze this in terms of confidence 
intervals. For the four observed events, with zero back- 
ground, the Poisson signal mean 99% confidence interval 
is 0.82 — 12.23 |44|. Thus our observed mean for random 
events, 0.81 ±0.01, falls at the lower edge of this interval, 
yielding a 1% probability for a chance occurrence. 

We now discuss the implications of our results. Clearly, 
spatial correlation analysis with a well defined and large 
sample of CR positions ought to provide the key to the 
identification of EHECR sources. The result embodied 
in Eq. is not compelling enough to definitively rule 
out chance probability as generating the correlation of 
the observed events with the candidate sources, but it is 
suggestive enough to deserve serious attention in anal- 
yses of future data. Besides, it should be stressed that 
the starburst hypothesis predicts a spectrum which is ap- 
proximately a smooth power law between 10 18 ' 7 eV and 
10 20 ' 5 eV, in very good agreement with that reported 
by the AGASA Collaboration In addition, we note 
that a medium mass nucleus fits the shower profile of the 
highest energy Fly's Eye event quite well [45j . Moreover, 
the high muon density observed in the Yakutsk event also 
favors a nucleus primary [25j | . 



LOW BUDGET ANISOTROPY TARGET FOR 
THE HIGHEST ENERGY COSMIC RAYS 

The superior angular and energy resolution of the 
Pierre Auger Observatory 0] will allow the high end 
of the energy spectrum and the CR arrival directions to 
be measured with unprecedented precision. The proto- 
cols for testing anisotropy detection claims are being re- 
stricted to hypotheses that must be specified a priori in 
order to ensure: (1) that the sample is not (inadvertently) 
devised to suit a special hypothesis after a preliminary 
study of the data, and (2) that the number of potential 
sources is not so large that the criterion for negating the 
null hypothesis for any source is reasonable |47j . Since 
the budgeting for candidate sources of anisotropy is very 
constrained (total random probability after accounting 
for all relevant trials is less than 0.001) particular em- 
phasis should be put - in our view - to those models 
in which there is no new physics involved, and a plau- 
sible astrophysical mechanism is suggested as the origin 
of some or all events. Some examples are: the above 
mentioned Centaurus A, nearby quasar remnants |48| . 
and luminous infrared galaxies |49j ■ In this direction, we 
believe the reasons we just listed above are sufficient to 
encourage the Auger community to search for evidence 
of the starburst model in forthcoming measurements. 

We now estimate the sensitivity of Auger to our model. 
The event rate for the Southern Auger Observatory (a de- 
tector with aperture A = 7000 km 2 sr above 10 19 eV, and 
angular resolution less than 1.5° with 8 < 20°), assum- 
ing extrapolation of AGASA flux {E 3 J(E) « 10 24 5 eV 2 



sr 



H3) up to 10 20 - 5 eV, is given by 
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where E x = 10 20 3 eV and E 2 = 10 20 5 eV. We now con- 
sider a 5-year sample of 25 events, and note that for the 
energy range under consideration the aperture of Auger 
is mostly receptive to cosmic rays from NGC253. We 
allow for different possibilities of the effective reduction 
of the cone size because of the Galactic magnetic field 
biasing discussed previously. In Fig. [3] we plot contours 
of constant probabilities (P = 10~ 4 , 10~ 5 ) in the two- 
dimensional parameter space of the size of the cone (as a 
fraction of the full 50° circle) and the minimum number of 
events originating within the resulting effective solid an- 
gle. Several important conclusions may be drawn. First, 
there is very little sensitivity of the results to the size 
of the cone, the variation is less than 20% for a 50% re- 
duction in the cone size. Secondly, the model predicts 
that after 5 years of operation, all of the 25 highest en- 
ergy events would be observed in the aperture described 
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FIG. 3: Curves of constant probabilities in the two- 
dimensional parameter space denned by the size of the cone 
and the minimum number of events originating within the 
resulting effective solid angle. 



above. From Fig. [3] we can see that even if 7 or 8 are ob- 
served, this is sufficient to rule out a random fluctuation 
at the 10 -5 level. Thus, the disproof of the starburst hy- 
pothesis can be achieved at a very small cost, < 1CP 5 out 
of a total 10 -3 to the Auger probability budget. Current 
preliminary assignments for other hypotheses are on the 
order of 10" 4 E3. 



CONCLUSION 

We have made a definite prediction for future obser- 
vations at the Auger Observatory: if the origin of CRs 
above 10 18 7 eV are nearby starburst galaxies, the in- 
coming CR flux will show a strong dipole anisotropy in 
the harmonic decomposition at energies beyond 10 20 3 eV. 
Because of its well-defined prediction, the model can 
be tested at the 5<r level in five years of running at 
Auger. Therefore, we strongly recommend that the 
Auger Collaboration take into account the next-door 
galaxy NGC253 in their first anisotropy prescription for 
super-GZK CRs. The confirmation of the starburst hy- 
pothesis would provide, as spinoff, direct evidence for the 
global structure of the Galactic magnetic field. 



Acknowledgments 

The work of L.A.A. and H.G. has been partially sup- 
ported by the US National Science Foundation (NSF), 
under grants No. PHY-0140407 and No. PHY-0073034, 
respectively. The work of D.F.T. was performed under 
the auspices of the U.S. Department of Energy by Uni- 
versity of California Lawrence Livermore National Labo- 
ratory under contract No. W-7405-Eng-48. 



[2. 



[9 
[10 



[11 
[12 

[13] 



K. Greisen, Phys. Rev. Lett. 16 (1966) 748; G. T. Zat- 
sepin and V. A. Kuzmin, JETP Lett. 4, 78 (1966) [Pisma 
Zh. Eksp. Teor. Fiz. 4, 114 (1966)]. 

T. Stanev, R. Engel, A. Mucke, R. J. Protheroe 
and J P. Rachen, Phys. Rev. D 62, 093005 
(2000) |arXiv:astrcTp h/0003484 ; F. W. Stecker and 
M. H. Salamon, Astrophys. J. 512, 521 (1992) 



arXiv:astro-ph/9808110 


; L. N. Epele and E. Roulet, 


JHEP 9810, 009 (1998) 


arXiv:astro-ph/9808104 




J. W. Elbert and P. Sommers, Astrophys. J. 441, 151 



(1995) arXiv:astro-ph/9410069 . 

The sphere within which a source has to lie in order to 
provide us with EHECRs if the primaries are subject to 
the GZK cutoff phenomenon. 

M. Ave, J. Knapp, J. Lloyd-Evans, M. Marchesini and 

A. A. Watson, arXiv:astro-ph/0112253 

N. N. Efimovi T. A. Egorov, A. V. Glushkov, 

M. I. Pravdin and I. E. Sleptsov, in Astrophysical Aspects 

of the Most Energetic Cosmic Rays, (Eds. M. Nagano and 

F. Takahara, World Scientific, 1991), p. 20. 

D. J. Bird et al. [HIRES Collaboration], Astrophys. J. 

424, 491 (1994). 

T. Abu-Zayyad et al. [High Resolution Fly's Eye Col- 
arXiv:astro-ph/0208243 T. Abu-Zayyad 



laboration] 
et al. [High 
|arXiv: astro-ph /020830 1 
J. N. Bahcall and E ~ 
(2003) |arXiv:hep-ph/0206217| . 
M. Takeda et al, Phys. Rev 
arXiv:astro-ph/9807193 ; 



Resolution Fly's Eye Collaboration] 



Waxman, Phys. Lett. B 556, 1 



[14] 



Lett. 81, 1163 (1998) 
N. Hayashida et al., As- 
trophys. J. 522, 225 T1999) |arXiv:astro-p h/0008102 ; 
M. Takeda et ffll., |arXiv:astro-ph/0209422| 
G. R. Farrar and T. Piran, |arXiv:astTo -ph/0010370 
M. Blanton, P. Blasi and A. V. Olinto, A stropart. Phys. 
15, 275 (2001) arXiv:astro-ph/0009466 . 
Y. Uchihori, M. Nagano, M. Takeda, M. Teshima, 
J. Lloyd-Evans and A. A. Watson, Astropart. Phys. 
13 (2000) 151 arXiv:astro-ph/9908193 ; P. G. Tinyakov 
and I. I. Tkachev, JETP Lett. 74, 1 (2001) [Pisma Zh . 
Eksp. Teor. Fiz. 74, 3 (2001)] [arXiv:ast ro-ph/0102101 ; 
L. A. Anchordoqui, H. Goldberg, S. Reucroft, 
G. E. Romero, J. Swain and D. F. Torres, Mod. Phys. 
Lett. A 16, 2033 (2001) arXiv:astro-ph/0106501 . 
E. Waxman, K. B. Fisher and T. Piran, Astrophys. J. 
483, 1 (1997) arXiv:astro-ph/9604005 . 



7 



[15] See e.g., P. Bhattacharjee and G. Sigl, Phys. Rept. 
327, 10 9 (2000) |arXiv:astro-p h/981 1011 ; S. Sarkar, 
arXiv:hep-ph/0202013 P. Bhattacharjee and G. Sigl. 
Prepared for Physics and Astrophysics of Ultrahigh- 
energy Cosmic Rays, Meudon, France, 26-29 Jun 2000; 
G. Sigl, arXiv:astro-ph/0210049 

[16] See e.g., R. J. Protheroe and T. Stanev, Phys. Rev. 
Lett. 77, 3708 (1996) [E rratum-ibid. 78, 3420 (1997)] 
arXiv:astro-ph/9605036 ; F. Csikor and Z. Fodor, light 
gluino window," Phy^ Rev. Lett. 78, 4335 (1997) 
arXiv:hep-ph/9611320 ; L. Anchordoqui, H. Goldberg, 
T. McCauley, T. Paul, S. Reucroft and J. Swain , 
Phys. Rev. D 63, 124009 (2001) arXiv:hep-ph/0011097 ; 
L. A. Anchordoqui, J. L. Feng, H. Goldberg and 
A. D. Shapere, Phys . Rev. D 66, 103002 (2002) 
arXiv:hep-ph/0207139 . 

[17] M. Lemoine, G. Sigl and P. Biermann, 
arXiv:astro-ph/9903124 G. R. Farrar and 

T. Piran, Phys. Re v. Lett. 84, 3527 (2000) 
[arXiv:astro-ph/9906431| . 

[18] L. Anchordoqui, H. Goldberg, S. Reucroft and J. Swain, 
Phys. Rev. D 64, 123004 (2001) |arXiv:hep-ph/0107287| . 

[19] K.-T. Kim, P. P. Kronberg, G. Giovannini and T. Ven- 
turi, Nature 341, 720 (1989). Fields of O(fiG) are also 
indicated in a more extensive study of 16 low redshift 
clusters. T. E. Clarke, P. P. Kronberg and H. Bohringer, 
Astrophys. J. 547, Llll (2001). 

[20] P. P. Kronberg, Rept. Prog. Phys. 57, 325 (1994); J. P. 
Vallee, Astrophys. J. 360, 1 (1990) 

[21] D. Ryu, H. Kang, and P. L. Biermann, Astron. Astro- 
phys. 335, 19 (1998). 

[22] The lower bound comes from the analysis in, L. A. An- 
chordo qui and H. Goldberg, P hys. Rev. D 65, 021302 
(2002) arXiv:hep-ph/0106217 . 

[23] P. Blasi and A. V. Olinto Phys. Rev. D 59 (1999) 023001 
arXiv:astro-ph/9806264 . 

[24] L. A. Anchordoqui, H. Goldberg and T. J. Weiler, Phys. 
Rev. Lett. 87, 081101 (2001) arXiv:astro-ph/0103043 . 

[25] N. N. Efimov, N. N. Efremov, A. V. Glushkov, 
I. T. Makarov and M. I. Pravdin, in Astrophysical As- 
pects of the Most Energetic Cosmic Rays, (Eds. M. 
Nagano and F. Takahara, World Scientific, 1991), p. 434; 
E. E. Antonov et al, JETP Lett. 69, 650 (1999) [Pisma 
Zh. Eksp. Teor. Fiz. 69, 614 (1999)]. 

[26] D. J. Bird et al, Astrophys. J. 441, 144 (1995). 

[27] N. Hayashida et al, Phys. Rev. Lett. 73, 3491 (1994). 

[28] N. Sakaki et al, in Proceedings of the 27th Interna- 
tional Cosmic Ray Conference ICRC2001 (Copernicus 
Gesellschaft, 2001), p.337. 

[29] C. Isola, M. Lemoine and G. Sigl, Phys. Rev. D 65, 
023004 (2002) |arXiv:astro-"p h/0104289 . 

[30] C. Isola and G. Sigl, arXiv:astro-ph/0203273 

[31] L. A. Anchordoqui, G. E. Romero and J. A. Combi, Phys. 



Rev. D 60, 103001 (1999) arXiv:astro-ph/9903145 . 

[32] G. Bertone, C. Isola, M. Lemoine and G. Sigl, 
arXiv:astro-ph/0209192 

[33] R. N. Manchester, Astrophys. J. 188, 637 (1974); R. C. 
Thomson and A. H. Nelson, Mont. Not. Roy. Astron. Soc. 
191, 863 (1980); J. L. Han and G. J. Qiao, Astron. Astro- 
phys. 288, 759 (1994); C. Indrani and A. A. Deshpande, 
New Astron. 4, 33 (1998). 

[34] R. Beck, A. Brandenburg, D. Moss, A. Shukurov, and D. 
Sokoloff, Ann. Rev. Astron. Astrophys. 34, 155 (1996). 

[35] J. L. Han, R. N. Manchester and G. J. Qiao, 
Mont. Not. Roy. Ast ron. Soc. 306, 371 (1999) 
[arXiv:astro-ph/990310l"1 . 

[36] J. P. Vallee, Astrophys. J. 366, 450 (1991). 

[37] T. Stanev, Astrophy s. J. 479, 290 (1997) 
[arXiv:astro-ph/9607086"l ; D. Harari, S. Moller- 
ach and E. RouletT JHEP 9908, 022 (1999) 
[arXiv:astr o-ph /990630 9"!. 

[38] R. Beck, arXiv:astro-ph/0012402 

[39] E. Waxman and J. Miralda-Escude, Astrophys. J. 472, 
L89 (1996) arXiv:astro-ph/9607059 . 

[40] J. Alvarez- Muniz, R . Engel and T. Stanev, 
arXiv:astro-ph/0112227 

[41] G. E. Romero, P. Benaglia and D. F. Torres, Astron. 
Astrophys. 348, 868 (1999) |arXiv:astro -ph/9904355 ; 
G. Sigl, D. F. Torres, L. A. A"nchordoqui and 
G. E. Romero, Phys Rev. D 63, 081302 (2001) 
[arXiv:as tro-ph /0008363] . 

[42] The standard deviation of the 4500 trial sample is 0.8, 
giving a standard deviation of 0.01 for the mean. 

[43] Because of constraints inherent in partitioning events 
among clusters, the distributions are very close 
to, but not precisely Poisson. H. Goldberg and 
T. J. Weiler, Phys. Rev. D 64, 056008 (2001) 
[arXiv:as tro-ph /0009378] . 

[44] G. J. Feldman a nd R. D. Cousi ns, Phys. Rev. D 57, 3873 
(1998) |arXiv:physics/9711021| . 

[45] F. Halzen, R. A. Vazquez, T. Stanev and H. P. Vankov, 
Astropart. Phys. 3, 151 (1995); L. A. Anchordo- 
qui, M. Kirasirova, T. P. McCauley, S. Reucroft 
and J. D. Swain, Phy s. Lett. B 492, 237 (2000) 
[arXiv:astro-ph/0007403"l . 

[46] See, for instance, L. Anchordoqui, T. Paul, S. Reucroft 
and J. Swain, jarXlydiep -ph/0206072 

[47] R. Clay, B. Dawson, B. Fick, M. Roberts, and P. Som- 
mers, GAP Note 2001-047. Electronically available at 
http : / /www . auger . org/admin/GAP_ Notes/index . html 

[48] D. F. Torres, E. Boldt, T. Hamilton and 
M. Loewenstein, Phys. Rev. D 66, 023001 (2002) 
[arXiv:astro-ph/0204419| . 

[49] A. Smialkowski, M". Giller and W. Micha lak, J. Phys. G 
28 (2002) 1359 |arXiv:astro-ph/0203337l . 



